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The edge states of monolayer and bilayer graphene ribbons both with/ without spin-orbit coupling 
and decoration of edges are studied by means of exact diagonal method. We show that the new 
extra edge states can be created by the decoration of edge layers. The fundamental picture of the 
competition among the Rashba, intrinsic spin-obit coupling (SOC), bias voltage on energy spectrum 
to effect on the energy spectrum of bilayer graphene and the phase of topological insulators is 
reported. 
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I. INTRODUCTION 

Recently, graphene^i^ has greatly attracted multifari- 
ous research activities owing to its rich and unexpected 
phenomena, also its remarkable promise for future nano- 
electronics application. It is an unique two dimensional 
(2D) system with one atomic layer in current physics. 
Its honeycomb structure provides some unconventional 
physics behaviors. Different from the usual parabola en- 
ergy spectrum of excitation in usual semiconductor and 
Fermi-liquid, the graphene has a Dirac-like linear, mass- 
less in low energy spectrum of excitation (up to ener- 
gies of the order of 1000 K) and the pseudospin degen- 
eracy that make the graphene to be a truly 2D "rela- 
tivistic" electronic system in solid and can be described 
by a low energy 2 + 1 dimensional effective massless 
Dirac theory^^ with a Fermi velocity vf ^ lO^m/s. Its 
massless linear energy spectrum presents unconventional 
physical behaviors such as its Integer Quantum Hall Con- 
ductance (IQHC) predicted theoretically^^^ and verified 
experiment all}J^i^, tunnelling through a potential bar- 
rier in graphene^ and others. Rather different from 2D 
Electron Gas (2DEG), graphene has two Dirac points 
in its Fermion spectrum which are at the Fermi energy 
for undoped system. In the presence of a perpendicu- 
lar magnetic field, the energy spectrum of graphene con- 
tains positive and negative Landau levels, as well as a 
zero energy Landau level. It is quite different from the 
standard 2DEG. Thus IQHC of graphene shows quite 
different behavior from conventional IQHC^^^^^. Ones 
have known that a wide class of topological phases is 
tightly related with the existence of edge states. IQHC 
and Quantum Spin Hall Conductance (QSHC) are two 
example^"—. Different from IQHC whose presence fun- 
damentally relies on the breaking of the Time-Reversal 
Symmetry (TRS) via applied external magnetic field or 
intrinsic SU(2) gauge flux—, TRS is retained in Quan- 
tum Spin Hall Effect (QSHE) that ensures the gapless of 
edge states and robust to non-magnetic disorder owing 
to Kramers degeneracy. It turns out that a Z2-valued 
topological invariant associates with QSHE^^j^i^ — . A 
theoretical model of graphene with SOC was suggested 
to explore the topological phase of QSHE by Kane and 
Mele^j^. In their work, following two kinds of SOC were 



introduced: 

i) . Intrinsic SOC provided from the next nearest neigh- 
bor (n.n.n.) spin-dependent hopping term that is some- 
thing similar to Haldane's picture^^. This kind of SOC 
can generate a bulk energy gap^^ and convert graphene 
from a 2D gapless semiconductor to a topological insu- 
lator that presents the QSHE^^^^^ and robust to non- 
magnetic disorder owing to TRS^-^^; 

ii) . Rashba SOG^^^i^"— induced by structure in- 
version asymmetry. It has been known that QSHE 
will be gradually destroyed as invasion of the Rashba 

However, according to the work by Y. Yao et al.^^, the 
SOC {10~^meV) in graphene is much smaller than the 
one suggested by Kane's modeU^. Therefore the SOC in 
graphene has been neglected in most works. However, 
the fundamental picture of QSHE in theoretical model 
with artificial SOC in graphene and clear understanding 
the competition among the Rashba, intrinsic SOC and 
other possible external controllable parameters to effect 
on the phase of QSHE is very importan t ^^i^^i^^ . 

In the physics, there might be some correspondence 
between bulk and edge states in topological states^ like 
the IQHC and QSHC, the study of edge states has be- 
came a hot topic recently. To study of the edge states of 
monolayer and bilayer graphene naturally becomes one 
of the research subjects. It has been found that a non- 
dispersive zero mode of the edge states (namely a fiat 
band in zero energy) appears in the graphene ribbons 
(GR) with zigzag edges (GR-ZZE). If the n.n.n. hop- 
ping t'{^ O.lt, here t is the energy of nearest neighbor 
(n.n.) hopping) is taken into account, its zero mode of 
the edge states will be changed from the non-dispersive 
to dispersive^i^. It has been also shown that no edge 
state exists in GR with armchair edges (GR-ACE). 

In bilayer graphene (BG), the spectrum of low energy 
excitation becomes massive, but it still keeps gapless^^. 
When an electric bias field is applied between two layers 
of BG, its gapless spectrum can be opened and the gap 
opened can be tuned by the bias voltage^"—. Such tran- 
sition from gapless to a gapful semiconductor has been 
observed in experiement^^. It was shown that the edge 
states do exist in BG-ZZE^^^, but not in BG-ACE. 

In terms of exact diagonal method, this paper would 
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give more detail studies on the edge states for monolayer 
and bilayer graphene ribbons (BGR) with/without SOC 
as well as the decoration of edges (DE). When the hop- 
ping constants connected to the edge sites are chosen as 
a different value from ones in bulk, it is found that some 
extra edge states can be created for both GR-ZZE and 
GR-ACE. It may suggest a way to manipulate the trans- 
port behavior of GR by means of decorating the edge 
layers. In our calculations, we show that different be- 
haviors even both the intralayer Rashba SOC and the 
interlayer Rashba go against the phase of topological in- 
sulators. The invasion of intralayer Rashba SOC in BG 
reduces the gap in bulk energy spectrum induced by in- 
trinsic SOC down to zerc^^, but the interlayer Rashba 
SOC will open a gap for the band of edge states fulfilled 
in the bulk gap region that drives such topological "edge 
semimetal" back to be an insulator. 

The present paper is organized as follows. In Sec. [Ill 
we summarize the results of edge states in graphene with 
and without DE, but no SOC is considered. The effect 
of SOC on the energy spectrum, particularly on the edge 
states, of GR is reported in Sec. IIIII Sec. IIVI presents 
our results on the edge states for BGR. Finally, we give 
a brief conclusion in sec. |Vl 

II. EDGE STATES IN GR WITHOUT SOC 

The electronic transport in graphene is mainly con- 
tributed hy Pz — 7T orbit als. A tight-binding approxima- 
tion (TBA) has been well applied and a theoretical model 
Hamiltonian is 

n,= Y^tijcicj+ J2 iii4ci + A^Ci4ci, (1) 

where the subscript in first summation represents 
a pair of n.n. lattices, <^i,j^ in second summation 
means a pair of the n.n.n. lattices, the hopping con- 
stants in bulk are tij{= t ~ 2.8eV) and t[j{= t' ^ O.lt). 
In our theoretical model, the n.n. hopping constants tij 
connected to edge sites are taken as different from 
to) to roughly approach some effective decoration 
of edges such as the softening of edge structure, atomic 
substitution and others. The last term in Hamiltonian 
describes the staggered potential (A 7^ 0,^i = ±1) added 
on two sublattices that can result the symmetry breaking 
in exchanging two sublattices A and B. Above Hamilto- 
nian can be diagonalized exactly for ribbons with differ- 
ent edge geometries such as bearded, zigzag and armchair 
(see Fig. [U^a)). Thus, the eigenvalues and eigenstates can 
be calculated numerically. Then the energy bands and 
the local density of states (LDOS) are obtained. In first, 
we repeat some results again which have been reported 
in the literatures —i^^i^ — for the comparison with our 
systems having some edge decoration. 

Comparing with the bulk energy spectrum of graphene, 
an extra non-dispersive zero mode of edge state appears 
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FIG. 1: Schematic illustration of lattice structure of mono- 
layer (a) and bilayer (b) GR. 



for zigzag-bearded edged graphene ribbon (GR-ZBE) 
(see Fig. [IJa)). When the n.n.n. hopping ^ 0.1 is 
added, the energy spectrum is shift down and the edge 
mode in GR-ZBE are changed from the non-dispersive 
(see Fig. Efb)) to dispersive (see Fig. Efd)). In contrast 
to GR-ZZE and GR-ZBE, no edge state can be found for 
GR-ACE (see Fig. EJa)) even n.n.n. hopping is included 
(see Fig. EJc)). The basic features in Fig. [2fa) and 
(c) are qualitatively the same where the linear dispersion 
near Dirac point K{K') in two cases are still preserved, 
but a small gap (see Fig. [2fa) and (c)) may appear due 
to the finite width of the ribbon. The metallic condition 
of GR-ACE, n = 3m + 2, is repeated, here n is chain 
number and m an integer. The GR-ACE should be an 
insulator when n ^ 3m + 2. 

The dispersion relation of the electron spectrum of 
both GR-ZZE and GR-ACE have been analytically stud- 
ied based on both continue Dirac equation^^^^^ and dis- 
crete model at TBA2ii^ii^-i& with proper boundary con- 
ditions. In both zigzag and bearded edged graphene, 
the zero mode of edge state appears and its wave func- 
tion vanishes at all sites of one sublattice and localize 
at another sublattice near the edges of the ribbon. For 
GR-ACE, no edge state is found. The dispersion rela- 
tions are shown in Fig. [2fa) and (b) that have been fre- 
quently appeared in literatures. For GR-ZZE, the edge 
states^-1^1^1^ exist only in region: 27r/3 < kx < 47r/3. 
However, the edge states exist only in — 27r/3 < kx < 
27r/3 for graphene ribbon with bearded-bearded edges 
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FIG. 2: (color online) The energy spectrum for graphene rib- 
bon with armchair edge (a) and (c), and zigzag-bearded edge 
(b) and (d). The hopping constants are t = 1.0,t^ = for 
(a)-(b), and t = 1.0, = 0.1 for (c) and (d). 



(GR-BBE). Hence, if the GR is zigzag-bearded edged^-, 
its edge states will spread to full region of kx G [— tt, tt] 
(see Fig. Efb) and (d)). 

When the A-B staggered potential A is applied, the A- 
B sublattice symmetry is violated and an bulk energy gap 
appears (see Fig. [3]). The magnitude of gap is approxi- 
mately proportional to the strength of staggered sublat- 
tice potential A. As we know, if the sublattice at one edge 
of GR-ZZE or GR-BBE is A(B), the one at another edge 
of GR must be B(A). But GR-ZBE is different, the type 
of atoms at two edges must be same. Thus, the energy 
spectrum of edge states shows something difference that 
have been shown in Fig. [3l 

From the physics, it is able to make some decoration 
to the atoms near edge surface^^, here we simply take an 
effective different hopping constant to in connecting the 
edge sites to replace the bulk value t. In the case, the 
spectrum of edge states clearly shows the changes. Some 
new extra edge states are created when to ^ t (see Fig. 
|4]and[5j). Fig. |4] gives the energy spectrum of GR-ACE 
with to 7^ where (a), (b), (c) corresponds to to = 0.2, 
to = 0.5 and to = 0.8, respectively. It is found that some 
new extra edge states are found at all kx points across 
the gap between conducting and valence bands. Chang- 
ing the ratio to/t from 0.2 to 1.0, the new extra edge 
states are changed dramatically. In the limit of to = t, 
the edge states fully disappear. It may suggest that the 
edge states can be manipulated by edge decoration. Very 
similar picture is also obtained for GR-ZZE, GR-BBE 
and GR-ZBE, some additional edge states also emerge in 
its energy spectrum (see Fig. [5j). The LDOS is attached 
in Fig. Ill^d) and Fig. ^d) where the LDOS show the 
exponent decay away from two edges of GR. Fig. [5fa) 
describes the energy spectrum of the edge states for GR- 
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FIG. 3: (color online) The energy spectrum of the graphene 
ribbon with different edges: (a) armchair- armchair, (b) 
zigzag-zigzag, (c) bearded-bearded, and (d) zigzag-bearded. 
All parameters are fixed at to —t= 1.0, t' = 0.1 and A = 0.1 
for (a)-(d). 



ZZE with edge decoration where the extra edge states 
can be found in the rest region of — 27r/3 < kx < 27r/3. 
Similar feature can be found for GR-BBE with edge dec- 
oration from Fig. Efb) and for GR-ZBE with same edge 
decoration from Fig. [5fc). The LDOS of edge states are 
plotted in Fig. Efd). All extra edge states disappear in 
the limit of to t. 



III. EDGE STATES IN GR WITH SOC 

Now we turn to study the tight-binding model of 
graphene with Rashba and intrinsic SOC that was first 
proposed by Kane and Mele^^^^^. The model Hamilto- 
nian of the graphene with SOC reads: 

^ iXn 4ai^ ^ dij)^Cj^, (2) 

where Hi is the free part of hopping terms which has 
been described in Eq. ([T]), the second term in Eq. ([2]) 
is the intrinsic SOC which is induced by spin-dependent 
n.n.n. hopping where Uij = {2/\^){di x d2)z = ±1, j) 
the index of sites, and (di, ^2) the unit vectors along the 
two n.n. bonds from site i to j, Pauli matrix in z direc- 
tion. Constant Xso is the intrinsic SOC and can be esti- 
mated in the range from ^ 0.001 to ~ O.OSmeV in terms 
of perturbation theory or the first-principles calculations. 
The third term in Eq. (|2]) is Rashba SOC induced by a 
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FIG. 4: (color online) The dispersion relation for graphene 
ribbon with armchair edge (a)-(c) and LDOS (d) correspond- 
ing to (b). The parameters: t = 1.0, = 0.1 for (a)-(c), 
to = 0.2 for (a), to = 0.5 for (b), and to = 0.8 for (c). 



FIG. 5: (color online) (a) is the energy spectrum of graphene 
ribbons for zigzag-zigzag edge , (b) for bearded-bearded edge, 
(c) zigzag-bearded edge and (d) for LDOS of the edge states 
in (c). The parameters are fixed at t = 1.0, t^ = 0.1, to = 0.5 
for (a)-(c). 



perpendicular electric field and n.n. hopping. The cou- 
pling constant A^^ can be experimentally determined by 
the spin angle-resolved photoemission spectroscope^^ or 
the measurement of spin relaxation^^^^^ in graphene. 



The effects of Rashba and intrinsic SOC on energy 
gap and edge states of the graphene with edge deco- 
ration are presented in Fig. [6l From Fig. [6l^a), one 
can see an energy gap and some extra edge states across 
the gap area where the spin degeneracy has been split- 
ted except the crossing points that is protected by TRS. 
Thus, it presents a topological insulator which the bulk 
is insulator, but becomes a metal at surface when an 
edge is inserted no mater what the geometric structure 
of edge is, armchair or zigzag edge anyway. Such degen- 
eracy resulted from TRS is robust against non-magnetic 
disorders^^^^^^^'^. If Rashba SOC Xr exists only, there is 
no energy gap opened so it is not able to be an insulator. 
Difference from Rashba SOC, the intrinsic SOC is favor- 
able to induce an energy gap around Dirac point^^. More 
lager Xso yields more big energy gap (see Fig. [6l^d)). 
When the edge is inserted into bulk, some edge states 
with spin chiral are created acrossing the gap that re- 
sults the system to be a metal since the edge states are 
extended along the edge surface and fullfill the energy 
gap. Rashba SOC will depress the gap that is unfavor- 
able to the topological phase^^. The Fig. [6l^d) shows the 
dependence of gap energy on Rashba SOC Xr. There is 
a critical value of above that the energy gap will be 
closed^^. 






-.-A^^ = 0.01 




-•--13^ = 0.03 
-T- = 0.08 



















FIG. 6: (color online) The energy spectrum of armchair edged 
graphene ribbon (a)-(c), and (d) shows energy gap Eg is func- 
tion of Xr at different values of Xso (d). The parameters 
Xso = 0.05 for (a) and (c), Xr = 0.05 for (b), and Xr = 0.2 
for (c). 



IV. EDGE STATES IN BGR WITH/ WITHOUT 
SOC 



We have mentioned that the Rashba SOC in mono- 
layer graphene contributes a factor to deduce the energy 
gap^^ although intrinsic SOC is favorable for opening an 
energy gap^^. But, it has been known that the gap is very 
small (zero) in graphene. However, it is reported that the 
energy gap can be opened in BG by an external electric 
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FIG. 7: (color online) The energy spectrum of bilayer 
graphene ribbons with armchair edge (a) and zigzag edge 
(b)-(d). The parameters t = 1, t± = 0.13 for (a)-(d), 
= 0.1, = 0.11 for (a), Vb = 0.1 for (a), (b) and (d), 
and Xso = 0.015 for (c)-(d). 

field^^^^. In this section, we would study the effect of 
SOC, bias voltage on the energy spectrum of BGR. 

A paper^^^ reported some results for unbiased and bi- 
ased bulk systems of bilayer graphene with both the in- 
tralayer Rashba and intrinsic SOC as well as an inter- 
layer Rashba SOC. As a result, under certain conditions, 
the energy spectrum of unbiased BG presents a Mex- 
ican hat-like energy dispersion. In the presence of an 
intralayer Rashba SOC only, another interesting result is 
that the K{K') point splits into four distinct ones that is 
contrastive to the case in monolayer graphene where the 
splitting also takes place, but the low-energy dispersion 
at these points remains identical. 

Here, we focus on the study of edge state for BGR (see 
Fig. Wp)' The Hamiltonian is given by 

<i,j> ,a(35 

+ "Y" (4,a,lCi,a,l - 4,a,2^j,^,2) + h.C. (3) 

i,a 

where 

^0=t Y 4a,S^j^,5 ^t' Y 4a,S^j^,5- 




1.74 2.44 3.14 3.84 4.54 ~ 1.74 2.44 3.14 3.84 4.54 



FIG. 8: (color online) The energy spectrum of zigzag edged 
bilayer graphene ribbon with SOC and bias voltage. The 
parameters t = 1, t± = 0.13 for (a)-(d), Vb = 0.1 for (a) and 
(b), Xr = 0.03 for (a) and (c), = 0.03 for (b) and (d), and 
Xso = 0.015 for (c) and (d). 



In above equations, we have used the additional sub- 
script (J to denote the spin polarization of electrons 
(cr =t,i) and 6 the layer number {S = 1,2). The sec- 
ond and third terms in H are the intrinsic and intralayer 
Rashba SOC. The fourth term describes n.n.n. interlayer 
Rashba SOC induced by a bias voltage Vb across two lay- 
ers of BG, the fifth and sixth terms describe the n.n. and 
n.n.n. hopping of interlayer^^, and the last term is the 
electric potential yielded by the bias voltage Vb- 

In Fig. [71 we give the energy spectrum of the electrons 
for BGR- ACE (a) and BGR-ZZE (b)-(d). In our precise 
calculation, it is shown that no edge state can be found 
in BGR- ACE even for biased case (see Fig. [Il^a)). As we 
have known that linear dispersive spectrum near Dirac 
point in unbiased BG is changed to be parabolic due to 
interlayer coupling, but no gap is opened. If a bias field 
is applied between two layers of BG, a gap will be open 
in bulk energy spectrun>^ (see Fig. [Tl^a) and (b)). The 
edge states are found for BGR-ZZE and shown in Fig. 
Wijo) where two degenerated modes of edge states in two 
decoupled GR-ZZE now are splitted and mixed due to 
the interlayer coupling^. If the n.n.n. hopping is taken 
into account, the electron- hole symmetry will be broken. 
In Fig. [71(c) and (d), the bulk energy gap is opened by 
intrinsic SOC, and edge states are found in BGR-ZZE. 
The difference of figure (b) and (c) is that the bulk gap 
in (b) is induced by bias voltage Vb only, but (c) by 
the intrinsic SOC only. They show very different feature 
of edge states. Figure (c) very clearly shows two Dirac 
cones where the up-down shift of two Dirac cones is from 
the interlayer coupling. The TRS at crossing points a 
and c in figure (c) is preserved, but the degeneracy at 
two crossing points b and d in figure (c) are accidental 
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and easily splitted by inter layer Rashba SOC A;^ that 
is shown in Fig. [Sfd). The difference of Fig. W[c) and 
Fig. [Tl^d) shows that further increasing the bias field will 
reduce the bulk energy gap opened by intrinsic SOC of 
BGR-ZZE. It is interested that bias voltage is favorable 
to open bulk energy gap for BGR-ZZE without intrin- 
sic SOC, but becomes unfavorable to keep a gap for the 
BGR-ZZE system with intrinsic SOC^^ 

In FiglHl we report some results on the competition 
among the Rashba, intrinsic SOC and bias field. We 
show that both the intralayer and inter layer Rashba SOC 
can induce the spin polarization and decrease the energy 
gap (see Fig. [8fa) and (b)), here the bulk energy gap 
has been opened by bias voltage Vb- The difference of 
(c) and (d) points out that the interlayer Rashba SOC 
A;^ changes the edge states to be gapful. But the Dirac 
points are still preserved by TRS. 

V. CONCLUSIONS 

In terms of exact diagonal method, we first report our 
studies on the edge states for the GR with edge deco- 
ration. When the edges of GR has no decoration, per- 
fect boundaries, the edge states exist only for GR-ZZE, 
GR-ZBE and GR-BBE, but not for GR-ACE, even for 
biased/unbiased BGR-ACE. However, some extra new 



edge states can be created for all GR with decorated 
edges including GR-ACE. In second, the effects of both 
the intralayer Rashba and intrinsic SOC on the energy 
spectrum are studied in detail for GR and BGR. It is 
shown that the Rashba SOC provides the unfavorable 
factor to the topological insulator through the deduction 
of the bulk energy gap. However, the intrinsic SOC drives 
the graphene to open a gap in the bulk spectrum^^, cre- 
ates the edge states full filled in the gap when an edge 
is inserted, generates Dirac cones protected by TRS that 
is favorable to form a topological insulator. At last, we 
also show that the bias voltage on BG, as same as Rashba 
SOC, can change the energy spectrum through the deduc- 
tion of the energy gap. But the interlayer Rashba SOC 
in BG destroy the topological phase of QSHE through 
opening the gap in the band of edge state fulfilled in bulk 
gap region that is also unfavorable to the phase of topo- 
logical insulator. Thus we conclude that these factors 
are all not favorable to the formation of the topological 
insulator except the intrinsic SOC. 
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